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Оцінювання прогнозів при моделюванні 
зсувних процесів Південного берега Криму 


Досліджено фактори та моделі зсувного процесу Південного берега Криму та наведено розрахунки 
оцінок параметрів цих моделей. Обрано критерії оцінювання моделей та оцінювання прогнозів, за якими 
виділяється модель серед інших. Модифіковано методику Бокса-Дженкінса відбору моделей та аналізу 
прогнозу. За допомогою модифікованої методики проведено оцінку обраних моделей та отримано прогноз. 


Вступ 


Прогнозування зсувних процесів Південного берега Криму, що викликають ката- 
строфічні наслідки, є основною, однією з найнеобхідніших, але при цьому найскладні- 
ших задач планування та аналізу. Наукові роботи (|1-3|, які присвячено аналізу зсувних 
процесів Південного берега Криму, спрямовані на статистичні спостереження, карту- 
вання, геологічні дослідження, описання процесів в цілому та експертну оцінку фахівця- 
геолога; слабкою ланкою цих досліджень є використання лише регресійного аналізу, 
відсутність системного підходу, обмежене використання сучасних інформаційних 
технологій та інтелектуальних систем прийняття рішень. Під системним підходом 
будемо розуміти системну узгодженість управління ресурсами, які мають бути вкла- 
дені в протизсувні заходи або у відновлення пошкоджених ділянок, та управління 
безпекою на всіх ділянках, що використовуються, не лише відповідно до мети, задач, 
ресурсів та очікуваних результатів, але також оперативності і результативності взаємо- 
дії в реальних умовах нештатних і катастрофічних ситуацій (4), 5. Задача оперативного 
прогнозування зсувних процесів Південного берега Криму - це задача управління 
екзогенними процесами, пов'язана з мінімізацією вимушених витрат на відновлення 
ушкоджених об'єктів. Ця задача є складною внаслідок недостатньої формалізованості 
та структурованості, багатокритеріальності, залежності процесу від поточної ситуації, 
високого рівня невизначеності та обмеження в часі на прийняття рішення. 

Метою даної роботи є моделювання та розробка довгострокових прогнозів 
зсувних процесів Південного берега Криму, для чого будуються математичні моделі 
зсуву, які найкращим чином описували б ретроспективні значення даних спостережень, 
модифікується методика Бокса- Дженкінса, за допомогою якої оцінюються моделі та 
прогнози. 


Аналіз факторів, які приводять до зсувів 


Розглянемо уважніше фактори, які є чинниками зсувних процесів на Південному 
березі Криму. Сукупність таких факторів можна поділити на три групи за ступенем їх 
змінюваності в часі: постійні, повільно та швидко змінювані ||, |3|. До першої групи 
можна віднести геологічні та геоморфологічні умови, що визначають генетичні особли- 
вості та інтенсивність появи зсувів. Фактори другої і третьої груп визначають тенденцію 
їх розвитку і діляться на дві підгрупи: незалежні і похідні. До незалежних відносять 
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тектонічний рух 1 клімат, до похідних - гідрогеологічні умови грунтового масиву. 
Атмосферні опади, температура навколишнього середовища, сейсмічна активність, 
сонячна активність, шторми, бурі та інші природні явища - це незалежні явища, ЯКІ 
швидко змінюються і діють на усталеність грунтів та сприяють появі напружених 1 
послаблених зон. 

При вивченні зсувів до початку 20-х років ХХ століття (В.Ф. Пчелінцев, Н.Ф. Пог- 
ребов, А.Ф. Ніфантов, І.Г. Глухов та ін.) підземні води вважалися головною причиною 
зсувів, а їх дренаж - основним заходом боротьби з ними. У той час було побудовано 
велику кількість підземних споруд, що виводили воду, але бажаного ефекту вони не 
дали. Таким чином, явної залежності між наявністю підземних вод на схилах і виник- 
ненням зсувів немає |2). 

Для підтвердження цього можна навести такі факти: 

- не на всіх схилах, які мають глинисті породи 1 підземні води, виникають зсуви; 

- в багатьох районах зі схожими геологічними і гідрогеологічними вимогами 
зсувонебезпечні схили перемежаються з усталеними; 

- є багато старих і стародавніх зсувів, зміщення яких припинилися за наявності 
тих же водоносних горизонтів, як і в період їх активності; 

- відомі «сухі» зсуви, які виникали за відсутності водоносних горизонтів, а також 
в суху погоду. 

Іноді вплив підземних вод на усталеність схилів великий, а іноді - зовсім незнач- 
ний. Хоч методи оцінки кількісного впливу є наближеними, ми маємо їх враховувати. 

Водоносні горизонти впливають на умови загальної рівноваги схилу в результаті 
створювання ними добавлених сил |З: 

- ваги вільної води, яка є у водоносному горизонті; 

- гідродинамічного тиску води, що рухається у водоносному горизонті. 

З часом з'явилася інша точка зору («Генсхема берегоукрепительньх сооружений», 
1949; И.Б. Корженевский и др., 1961), за якою головна роль відводилася абразії, ерозії, 
перевантаженню верхньої частини схила обвалами та осипами. Пізніше експеримента- 
ми зі штучної активізації показано, що найбільш дійовим фактором виявилася підрізка 
схилу в низовій частині незалежно від її генезису (природна або техногенна). 

Кількість техногенних зсувів, що з'явилися в результаті підрізки схилу, більше 
зустрічається на природних схилах із крутизною від 100 до 150 (46,997) ії від 150 до 
200 (31, Т2)., Менша кількість зсувів (від 2,690 до 10,290) формується як на пологих 
схилах від 5 до 10, так і на крутих - від 20 - 25 до 25 -35. , що відображає особ- 
ливості не тільки граничної рівноваги схилів, але Й їхнього освоєння (чим крутіше 
схил, тим рідше він використовується для розміщення об'єктів) |З . 

Вплив землетрусів на зсувні схили як чинник їх розвитку багатьма авторами 
відкидався. Останніми дослідженнями (безпосередньо для південноузбережного схилу) 
показано достовірний зв'язок сейсмічної і зсувної активності регіону впродовж року. 
Вплив сейсмічних чинників за певних сукупних умов (надмірне зволоження) може 
набувати не тільки допоміжного, але й першорядного значення. Тому, якщо не врахо- 
вувати цей чинник при аналізі режиму зсувної активності на ПБК, можна прийти до 
помилкових результатів, особливо при розгляді динаміки окремих зсувів. 

Будемо розглядати атмосферні опади як один з вирішальних факторів, що поро- 
джує зсуви на ПБК і визначає баланс вологи грунтового масиву. Велика кількість 
зсувів відбувається в період весняно-літнього періоду, коли в горах тане сніг або в 
період літніх злив, які відбуваються після тривалої засухи. Однак коефіцієнт кореляції 
становить достатню величину для початку моделювання лише при врахуванні опадів 
за два гідрогеологічні роки: попередній та поточний. 
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Головними факторами виникнення зсувів є такі види атмосферних опадів |31: 

- велика кількість одноразових опадів, тобто зливи або різке танення снігу; 

- невелика кількість опадів протягом тривалого часу (дощ, що мжичить від 3 до 
10 днів 1 більше). 


Методика прогнозування зсувних процесів 
Південного берега Криму 


У даній роботі запропоновано таку методику побудови моделей та побудови прог- 
нозів: всі можливі моделі із певного класу моделей аналізувати в три етапи |бЇ. 

На першому етапі аналізується якість моделі за такою схемою: 

1) коефіцієнт множинної детермінації: 


де у; - значення показника для /-./,..., п, у, - прогнозна величина, у - середнє; 
2) середній квадрат похибки: 
1 п 
2 СУД 
е 22У0, зб 
Па 
3) перевірка адекватності моделі за допомогою Е-критерію Фішера: 
п КУ з 
(пото 1) (5, - У) 
Е пе ізі . 
п ад 
тУ (у, - У) 
ізі 
4) перевірка на статистичну значимість моделі та її параметрів (і-критерій Стью- 
дента); 


5) перевірка на мультиколінеарність; 
6) перевірка на автокореляцію збурень (критерій Дарбіна- Уотсона), 


З (е, - еа) 


де р - коефіцієнт кореляції похибок моделі між собою; 

7) перевірка моделей на гетероскедастичність (тест Голдфелда та Квондта). 

За перерахованими критеріями обираються кращі моделі-кандидати та за їхньою 
допомогою будуються прогнози. 

На другому етапі оцінюється якість прогнозу за допомогою таких статистик: 

1) середня абсолютна похибка МАН (теап аб5оїшіе еттогу): 


І бе 
МАЕ дій 7 У; , 


ізі 


де у; - значення періоду, у; - прогнозна величина для 1 - 1,2,..., п; 
2) середня абсолютна відсоткова похибка МАРЕ (теап аб5оїиіе регсепіаєе егтог): 
мАрЕ З УРТІЇ лоді; 
Па і 
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3) середньоквадратична похибка КМ5Е (тоої теап 5дцагеа еттог): 


П - з 
ВММ5- - (у, - У) 
із! 


4) коефіцієнт нерівності Тейла : 


1 РОЮ; 
пн п 2075) | 
Пот айс ЗРІро об 
оо «ба 


5) для абсолютної 1 відносної абсолютної похибок обчислюються кількість у 
відсотках спостережень, які попали в інтервали |0; 0,50), Гс; оо), Г090; 109), 1509; со), 
а також їх максимальне значення для кожної моделі. 

На останньому, третьому, етапі проводиться оцінювання стійкості прогнозу, 
тобто значення всіх факторів, крім одного, фіксуються, а обраний фактор змінюється 
на Є середнє квадратичне відхилення (при цьому розраховується змінювання резуль- 
туючого показника у відсотках до його середнього значення). 


Оцінювання моделей і прогнозу зсувних процесів 
Південного берега Криму 


Таким чином, на першому етапі було відібрано серед інших п'ять моделей 
зсувних процесів Південного берега Криму, які будемо розглядати й надалі |71: 

- регресійна однофакторна модель (РО): АЗ(ї) - ад ч- аг Оці - 1); 

- авторегресійна модель (АР): АЗ(і) - ад ч аГ-АЗ(І - 1); 

- комплексна множинна лагова авторегресійна модель (ЛАР): 

АЗ(І) - а б агО( - 1) - а, СеА( - 1) 5 аз СоА(і - 1) б аг АЗ(І - 1) є ах (і - 1), 

де О - опади за два роки, СеД - сейсмічна активність, СоА - сонячна активність, 
АЗ - кількість активних зсувів за рік; 

- метод групового врахування аргументів (МГУА): 


М М М 
14 зач Ух); У ауд, -У У. не ЗР 
і-і уч іх) і-ї усі Кк) 


- метод комплексування аналогів (КА), за яким розраховується відстань між фак- 
торами історичних даних і факторами нового спостереження: 


ке У уза) -»тіпуіз1К, / 214, 
ра 


нехай Е,; - найменша відстань між факторами, Е,2 - наступна за нею, уд - величина у, 
що відповідає факторам з відстанню Е,/, а у,о - З Е;», тоді результат визначаємо за 
формулою: 


- Уа Ві о й ; 
вач» 


а! 


у 


На другому етапі проведемо оцінювання і порівняння обраних п'яти моделей 
за визначеними критеріями якості моделі і якості прогнозу (табл. 1). Наведена таблиця 
свідчить, що за якістю моделей та за якістю прогнозу модель МГУА - найкраща, за 
нею слідує ЛАР (найгіршою є регресійна однофакторна модель відносно опадів (РО) 
та модель комплексування аналогів (КА )). Проте жодну з наведених моделей відкида- 
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ти не варто тому, що вони по-різному реагують на різке збурення будь-якого з пара- 
метрів або всіх одразу, і тільки їх узагальнення дає найбільш наближений прогноз. 
Наприклад, модель комплексування аналогів в 4090 дає відносну похибку прогнозу, 


яка не перевищує 1097, і майже в 14905 спостережень перевищує 5097. 


Таблиця 1 - Аналіз якості моделі та якості прогнозу 


Якість прогнозу 
Моделі Абсолютна похибка Відносна похибка 
« 0,5с »с тах «1090 25090 тах 
АРИ1) 51,296 16,39 101,7 30,29 11,69 194,39 
РО(0) 55,89 23,39 120,9 39,59 18,69 210,69 
ЛАР 51,29 4,796 89,8 32,69 11,696 147,296 
КА 56,86 18,29 134 40,99 13,696 180,59 
МГУА 81,490 2,390 81,7 44,29 7,096 141,190 
Узагальнена 65,190 9,39 89,0 41,99 9,396 159,09 


Перевіримо таким чином кожну модель, додамо до цього критерію ще кількість 
абсолютних похибок, які належать указаним проміжкам, та максимальне значення 
абсолютної та відносної похибки (табл. 2). 


Таблиця 2 - Відсоток помилок, який належить указаному інтервалу 


Моне Якість моделі Якість прогнозу 
К2 Уе2/п | рМ/ СКВ | САВ | САВП | К-т Тейла 0 
АР(1) 0,501 1625 | 2,020 | 40,32 | 31,39 | 28,390 0,139 
РО(0) 0,422 1926 | 0,937 | 43,88 | 32,60 | 34,89 0,151 
ЛАР 0,643 1190 1,341 | 34,50 | 28,37 | 27,99 0,118 
КА 0,502 2206 1,483 | 46,39 | 32,19 | 27,29 0,167 
МГУА 0,805 648 2,301 | 25,46 | 18,98 | 18,99 0,087 
Узагальнена | 0,544 1093 1,500 | 42,19 | 24.36 | 24,19 0,135 
Наведена таблиця свідчить, що за якістю прогнозу за обома видами похибок 
модель МГУА - найкраща: 8190 від загальної кількості залишків становить менше за 


половину стандартного відхилення с спостережуваного показника, лише 290 -- переви- 
щує сб, максимальна абсолютна похибка становить 81,7, а відносна - 14190 (в той час, 
як для РО цей показник - 120,9 1 210,69 відповідно), 44790 від загальної кількості за- 
лишків становить менше 1090 відносної похибки і лише 795 перевищують 5090; за нею 
слідує ЛАР (найгіршою є регресійна однофакторна модель відносно опадів (РО) та 
модель комплексування аналогів (КА)). 

На третьому етапі проведемо оцінювання стійкості прогнозу, тобто значення 
всіх факторів, крім одного, фіксуються (береться середнє значення), а обраний фактор 
змінюється на - середнє квадратичне відхилення (при цьому розраховується зміна 
результуючого показника у відсотках до його середнього значення). Для обчислень 
використаємо програмний продукт РІ.55СС «Когесазипє ої І апдзде бог Зоціпегп Соазі 
ої Сгітеа» |8|, інтерфейс якого наведено на рис. І. 

Введемо спочатку середнє значення всіх факторів і обчислимо кількість активних 
зсувів у відсотках до середньої кількості. Додамо до всіх факторів середнє квадра- 
тичне відхилення та віднімемо його, результати обчислень помістимо в табл. 3. 
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Таблиця 3 - Середнє значення моделі для всіх факторів «с за 2010 рік 


РА АР ЛАР КА МГУА | Разом 
Хсер і 100,09 | 100,796 | 139,79о | 94,99 | 145,695 | 107,490 
Хсер ізо; 127,996 | 130,99 | 165,40 | 144,99 | 164,75 | 139,790 
Хсер і-ої 72,80 | 69,990 | 103,20 | 32,96 | 114,96 | 72,890 


Останній стовпчик показує узагальнене значення, яке розраховується як середнє 
зважене всіх отриманих прогнозів, вагами є коефіцієнти детермінації, поділені на їх суму. 


У 


2 2 2 2 2 
і Ура" Кра Я Удро Кар Удар Кар Ук ка В УурудКууд 


заг 


2 2 2 2 2 
Кр т Кр т Кр Р Ка Я Кигул 


4 Долгосрочньйй прогноз 
Вьіберите метод прогнозирования 
РА - регрессионньй анализ 
АР - авторегрессия 
ЛАР - лаговая авторегрессия 


МГУА - метод группового учета 
аргументов 
о Все методьї одновременно 


Введите даннье за прошльйй год: 


Сейсмическая активность 
Солнечная активность 


Количество активньхх оползней 


КА - комплексирование аналогов 


Вьполнить расчет 

| Назад | | Даннье | 

Результат показать 
в таблице 


на графике 


2010 


Количество осадков за два предьщущие года 1253 


4 


70 


135 


Вернуться и вьбрать другой тип прогноза 


1 Долгосрочньій прогноз 
Вьберите метод прогнозирования 


РА - регрессионньй анализ 

АР - авторегрессия 

ЛАР - лаговая авторегрессия 

КА - комплексирование аналогов 


МГУА - метод группового учета 
аргументов 


о Все методь одновременно 


мім  |МАХ |РА | (АР 


|ЛАР КА 


Назад | | Даннье | 
Результат показать 


в таблице 
на графике 


МГУА | 


Оползнізб 


264 | (136 | (137 
0,422 (0,501 


182 129 1192 
0,643 0,502 0,805 


4388 


Ожидаєемое количество оползней в следующем году (с учетом 
всех прогнозов) - 160 


Вернуться и вьібрать другой тип прогноза 


Рисунок І - Вікно вводу даних 1 результат прогнозування 


Табл. 3 свідчить, що найкращі моде; 


лі (МГУА 1 ЛАР), відмічені на попередніх 


етапах аналізу, поводять себе не зовсім коректно, тобто при зменшенні значень факторів 
вони показують збільшення досліджуваного показника. Наприклад, модель МГУА при 
збільшенні всіх факторів на середнє квадратичне відхилення дає прогноз, що перевищує 
середнє значенні майже на 6590, але при такому ж зниженні - все одно збільшується 
на 1496 (ЛАР - на 3,290). Це пов'язано з тим, що в ці моделі входить як фактор рік, на 
який робиться прогноз, тобто чим більший рік, тим більше результат. 

Проведемо порівняльний аналіз оцінювання прогнозу на 2010 рік двох найкращих 
моделей - МГУА 1 ЛАР. Підставимо всі середні значення для факторів, потім змінимо 
їх всі на єс, а потім лише один з факторів будемо змінювати, а інші дорівнювати- 
муть середньому. Результати розрахунків помістимо в табл. 4. 


Таблиця 4 - Оцінювання прогнозів на 2010 рік за двома моделями - ЛАР і МГУА 


2010 ЛАР | МГУА 2010 ЛАР | МГУА 
О(Хсері- 81) 14790 | 1949 О(Хсері1-51) 10496 13290 
Сед(Хсер2--52) 12190 | 15290 Сед(Хсер2-52) 10196 14090 
СоА(Хсер3--53) 11390 | 13790 СоА(Хсер3-53) 10896 16196 
АЗ(Хсера--54) 13890 | 15390 АЗ(Хсера-54) 10396 13390 
Середнє (У) 1309 | 15690 Середнє (У) 10490 1419 
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Наведена таблиця свідчить, що результуючий показник для моделі МГУА в серед- 
ньому збільшується на 5697 (від 3790 до 349), а замість зменшення в протилежному 
випадку зростає в середньому на 4190 (від 590 до 6190) залежно від обраного фактора. 
Модель ЛАР теж збільшується, але не так критично. 

За оцінками якості моделі та якості прогнозу модель МГУА є не тільки найкращою, 
а й випереджає інші на цілий порядок. Проте в табл. 3 1 4 доведено, що обидві моделі не 
є усталеними, вони мають великі відхилення при прогнозі на декілька кроків вперед. 

Таким чином, можна зробити висновок, що, дійсно, всі запропоновані моделі є 
адекватними до історичних даних 1 їх слід використовувати для прогнозування зсувних 
процесів Південного берега Криму, тобто не варто відхиляти одну чи декілька з них, 
а останній узагальнюючий результат можна розглядати як шосту модель. У табл. 112 
наведено оцінки цієї моделі. 

У табл. 5 наведено прогноз на 2005 рік, видно, що до фактичного значення най- 
ближчим є саме узагальнююче значення результуючого показника. 


Таблиця 5 - Прогноз на 2005 рік 


2005 РА АР ЛАР КА | МГУА | Узагальн. Фактичне 
прогноз 174 186 227 179 201 196 191 


Таким чином, наведені моделі використаємо для наглядного зображення прогнозу 
за допомогою пакета статистичної обробки інформації ЕМтіему8 (рис. 2): на чотирьох 
графіках (для моделей РА, АР, МГУА, ЛАР) дві верхні криві показують фактичне і 
прогнозне значення, а нижня - помилки прогнозування. 

я РУієму5 
Еїе ЕЕ Обіесів Мем» Ргосз Оціск Орбіоп5 Міпдом Неїр 


ни Едиаїіоп; ЕОРО Могікбіїе: ЦМТІТЕР Г. ГЕО) ла Бдиакіоп: ЕОРОС М/огікбїв: ПІМТІПЕР І- ех 


Мівмм| Ргосє | Обіесіз| | Ргіпі | Мате | Ртееге| Езіїгпаїе| Рогесазі| 5іаїз | Незійз Мівми | Рос Обіесіз| Ріїпі | Мате | Ртееге| Езіїгпаїе| Рогесазі | 5іаїз Весічіз | 


Незідизі Апизі Еінед Невідцаї Асіша! Еїнед 


на Едиаїіоп: ЕОМСЦА Могікбіів: ОМТІП.ЕЮ 1- ОЇХ| ма вдиагіоп: ЕОГАВ  М/огібїе: ОМТІПЕЮ 22) 


Мівми! Ргосе | Обіесіз | Риіпі | Мате! Ргееге| Езіїглаїе| Рогесагі | 5іаіз | Везідз Міеми! Ргосе | ОБівсіз| | Риїпі | Мате! Ргевге| Езіїглаїе| Рогесазі | 5іаїз | Везідз: 


Б5. 7ПО 75. 8п | 85. 9п 95 пп 


Рисунок 2 - Прогнозування зсувних процесів на Південному березі Криму 
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Отже, можна помітити, що помилки сягають особливо великих значень останніми 
роками (1993 - 2003). Це пов'язано більшою мірою з тим, що при Радянському 
Союзі після значних зсувів 1968 - 1969 років проводилися обстеження та укріплення 
зсувонебезпечних ділянок, після розвалу СРСР впроваджені заходи утримували схили 
від катастрофічних руйнацій, проте впродовж останнього часу на подібні роботи виді- 
ляється занадто мало коштів, деякі установи було скорочено або закрито, роботи прово- 
дилися більшою мірою з відновлення будов, а не попереджувального характеру. Все це і 
спровокувало сплеск активізації екзогенної активності. 

Таким чином, за отриманими моделями однофакторної регресії, авторегресії, ла- 
гової авторегресії, методу комплексування аналогів та методу групового урахування 
аргументів ми можемо розраховувати зсувну активність. Такий вид прогнозу доцільно 
використовувати при стратегічному плануванні, що надає розрахунки відповідних 
показників на наступні 1 - 5 років. 


Висновки 


У статті розглянуто моделювання зсувних процесів на Південному березі Криму 
на основі регресійного аналізу, методу комплексування аналогів і методу групового 
урахування аргументів, що дає можливість за отриманими моделями розраховувати 
зсувну активність. Проведено порівняльний аналіз цих моделей, а також на їх основі 
побудовано узагальнюючу модель. Таке оцінювання прогнозу доцільно використовувати 
при плануванні відповідних показників на наступні І - 5 років. 

Розроблено методику відбору і оцінювання моделей зсувних процесів Південного 
берега Криму, яка базується на стандартних методах перевірки якості регресійної 
моделі з урахуванням її адекватності історичним даним. На основі наведеної методики 
визначено, що обрані моделі коректно описують ці процеси. Результати прогнозування 
наведено у графічному вигляді. 
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Оцениваниєе прогнозов при моделированиий оползневьїх процессов Южного берега Крьшма 
Исследовань факторьг и модели оползневого процесса Южного берега Крьтма и приведеньт расчетьт 
оценок параметров зтих моделей. Вьбрань критерий оценивания моделей и оценивания прогнозов, 
по которьм вьщделяєтся модель среди других. Модифицирована методика Бокса-Дженкинса вьбора 
моделей и анализа прогноза. При помощи модифицированной методики проведена оценка вьтбранньх 
моделей и получен прогноз. 


У.М. Тагап 
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Чекегтіпайоп5 ої е5е глодеіз аге дїуеп Пеге. Стіегіа ої еуаїшайоп ої плодсіз апа Гогесазі5 Гог деїегтіпіпе 
Фе тодеі атопе офег тодеіз Ба8 Бееп 5еЇссіса. Те Вох-ІепКіп5 теїродоїову бог плодеі 58еЇесіїоп Гогесазі 
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